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1 Introduction 

Although the study of electron spin dynamics using pulse electron spin reso- 
nance is an established research field [T] , many theoretical questions regard- 
ing the microscopic mechanisms for reversible and irreversible decay of spin 
coherence remain open. Recently, the quest towards scalable quantum com- 
putation using electron spins [2] gave new impetus to pulse spin resonance, 
and sparked major experimental progress towards control and detection of 
individual electron spins in the solid state environment [3l |4j [6l [7] . 

The microscopic understanding of the mechanisms leading to electron spin 
energy and phase relaxation, and the question of how to control these pro- 
cesses is central to the research effort in spin-based quantum computation. 
The goal of theory is to achieve microscopic understanding so that spin coher- 
ence can be controlled either from a materials perspective (i.e., choosing the 
best nanostructure for spin manipulation and dynamics) or from the design of 
efficient pulse sequences that reach substantial coherence enhancement with- 
out a high overhead in number of pulses and energy deposition. (The latter is 
particularly important in the context of low temperature experiments where 
undesired heating from the microwave excitations must be avoided). 

One has to be careful in order to distinguish the time scales characterizing 
electron spin coherence. It is customary to introduce three time scales, T±, 
T 2 , and T 2 . For localized electron spins, these time scales usually differ by 
many orders of magnitude, because each is dominated by a different phys- 
ical process. T\ is the 1/e decay time for the spin magnetization along the 
external magnetic field direction. As an example, T\ of a phosphorus donor 
impurity in silicon is of the order of a thousand seconds at low temperatures 
and moderate magnetic fields (T = 4 K and B = 0.3 T) [8]. These long 
T\ 's are explained by noting that the spin-orbit interaction produces a small 
admixture of spin up/down states; the electron-phonon interaction couples 
these admixtures leading to tjt oc B 5 at low temperatures [HI [TO] ■ T\ is gen- 
erally long because a spin-flip in a magnetic field requires energy exchange 
with the lattice via phonon emission. The time scales T 2 * and T 2 are instead 
related to phase relaxation, and hence do not require transmission of energy 
to the lattice. Here T 2 is the 1/e decay time of the precessing magnetization 
in a free induction decay (FID) experiment (7r/2 — t — measure, where ir/2 
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denotes a spin rotation around the x axis). Hence T 2 * is the decay time of 
the total in-plane magnetization of an ensemble of spins separated in space 
or time (e.g., a group of impurities separated in space, or the time-averaged 
magnetization of a single spin, as discussed in section 12*751 below - ) . For a phos- 
phorus impurity in natural silicon, T 2 * « 20 ns due to the distribution of 
frozen hyperfine fields, that are time-independent within the measurement 
window of the experiment. The T 2 decay is reversible, because the ensemble 
in-plane magnetization is almost completely recovered by applying a spin- 
echo pulse sequence. In this review wc define as the 1/e decay time of a 
Hahn echo (ir/2 — t — it — t — echo). The irreversible decoherence time T2 is 
caused by uncontrolled time dependent fluctuations within each time interval 
r. For a phosphorus impurity in natural silicon, we have T2 ~ 0.3 ms 
four orders of magnitude longer than T 2 . 

The discussion above clearly indicates that the resulting coherence times 
are critically dependent on the particular pulse sequence chosen to probe spin 
dynamics. In section[2]we show that spin coherence can be directly related to 
the spectrum of electron spin phase fluctuations and a filter function appro- 
priate for the particular pulse sequence. The phase of a precessing electron 
spin is a sensitive probe of magnetic fluctuations. This gives us the opportu- 
nity to turn the problem around and view pulse electron spin resonance as a 
powerful tool enabling the study of low frequency magnetic fluctuations aris- 
ing from complex many-body spin dynamics in the environment surrounding 
the electron spin. 

A particularly strong source of magnetic noise arises due to the presence of 
nuclear spins in the sample. It is in fact no surprise that the dominant mecha- 
nism for nuclear spin echo [12] and electron spin echo decay [13l [11] has long 
been related to the presence of non-resonant nuclear spin species fluctuat- 
ing nearby the resonant spin. Nevertheless, the theoretical understanding of 
these experiments was traditionally centered at phenomenological approaches 
[TH HT5] , whereby the electron phase is described as a Markovian stochastic 
process with free parameters that can be fitted to experiment (this type of 
process has been traditionally denoted spectral diffusion, since the spin res- 
onance frequency fluctuates along the resonance spectrum in a similar way 
that a Brownian particle diffuses in real space). 

Recently, we embarked on an effort aimed at understanding the mech- 
anism of electron spin coherence due to nuclear spins from a fully micro- 
scopic point of view. In Reference [16] we developed a semiclassical model for 
electron spin echo decay based on the assumption that the relevant nuclear 
spin dynamics results from pair "flip-flops" , where the spin of two nuclei lo- 
cated close to each other is exchanged due to their mutual dipolar interaction 
(Fig.[T]). The flip-flop processes lead to fluctuations in the nuclear spin hyper- 
fine field seen by the localized electron (e.g., a donor impurity or a quantum 
dot in a semiconductor). The semiclassical theory is based on the assump- 
tion that each flip-flop can be described by a random telegraph noise process 
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Fig. 1. The electron spin of a donor impurity in silicon is sensitive to the magnetic 
noise produced by nuclear spins within its wave function. When two 29 Si isotopes 
are close to each other, their nuclear spin states may flip-flop due to their mutual 
dipolar interaction. These flip flop events produce time dependent fluctuations in 
the electron's hyperfine field, leading to phase relaxation and spin echo decay. 

(a phenomenological assumption), but with relaxation parameters that can 
be derived theoretically from a microscopic theory based on the nuclear spin 
dipolar evolution. Therefore this theory describes the irreversible decay of the 
effective hyperfine field produced by a pair of nuclear spins on the electron 
spin. Comparison with experiment [ITJ \17\ [lSl EH [2Q] suggested reasonable 
order of magnitude agreement for the 1/e echo decay time (within a factor 
of 3) but poor qualitative agreement for the time dependence of the echo 
envelope. The next step was to develop a full quantum theory for the nuclear 
spin dynamics affecting the electron spin. In reference [21] a cluster expan- 
sion method was developed to calculate echo decay due to the closed-system 
dynamics of a group of dipolar coupled nuclear spins, without any stochastic 
assumption about the nuclear spin dynamics. At lowest order in this cluster 
expansion the qualitative and quantitative agreement with experimental data 
was quite good. 

In section |4] we develop a fully microscopic theory for the nuclear spin 
noise spectrum arising due to pair flip-flops induced by the inter-nuclear 
dipolar interaction. This allows us to give an elegant and simple derivation 
of the lowest order cluster expansion results of Ref. [3T] and to interpret 
these results from the point of view of non-equilibrium statistical mechanics. 
The full noise spectrum is expressed as a sum of delta-function contribu- 
tions corresponding to isolated pair flip-flop transitions. We then show that 
irreversibility can be incorporated into the pair flip-flop processes by adding 
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broadening to these sharp transitions, in a mean-field like approach. Using 
the method of moments we are able to calculate these broadenings exactly 
(at infinite temperature). We show explicit numerical results for the noise 
spectrum affecting a donor impurity in silicon and compare the improved 
theory with echo decay experiments in natural |18j and nuclear-spin-enriched 
samples [19]. 

2 Noise, relaxation, and decoherence 

When the coupling between the spin qubit and the environment is weak, we 
may write a linearized effective Hamiltonian of the form 



where B || z is a static (time independent) magnetic field, 7 is a gyromagnctic 
ratio in units of (sG)^ 1 (we set Ti = 1 so that energy has units of frequency), 
a q for q = x,y, z are the Pauli matrices describing qubit observables and 
f] q {t) represents the environmental (bath) degrees of freedom. 

2.1 The Bloch-Wangsness-Redfield master equation 

In order to describe the long-time dynamics we may take the limit t — ► 
00. Such an approximation is appropriate provided t r c , where t c is a 
typical correlation time for bath fluctuations (later we will define r c properly 
and relax the long time approximation). In this case spin dynamics can be 
described by the Bloch-Wangsness-Redfield theory. The average values of the 
Pauli operator satisfy a Master equation (for a derivation see section 5.11 of 
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In Eqs. ([3]), (@|, (0 we assume (?)g(i)%'(0)) = for q ^ q' . These equations 
are the generalization of Fermi's golden rule for coherent evolution. From 
Eq. © we may show that the coherence amplitude |(a+)| = \{c x + ia y )\/2 
in a FID experiment decays exponentially with a rate given by 



In contrast, T\ is the time scale for (a z ) to approach equilibrium, i.e., T\ is 
the energy relaxation time. According to Eq. ([7]) we have T 2 * < 2T\. Note 
that Ti depends on the noise spectrum only at frequencies +7 2? and —jB, 
a statement of energy conservation. Positive frequency noise can be inter- 
preted as processes where the qubit decays from f to J. and the environment 
absorbs an energy quantum jB, while negative frequency noise refers to qubit 
excitation (from J, to j ) when the environment emits a quantum jB. The cor- 
relation time r c can be loosely defined as the inverse cut-off for S q (u>), i.e., 
for lu ^> 1/t c we may approximate S q (u>) ~ 0. 

The Master equation [Eq. (O] leads to a simple exponential time depen- 
dence for all qubit observables. Actually this is not true in many cases of 
interest, including the case of a phosphorus impurity in silicon where this ap- 
proximation fails completely (for Si:P the observed free induction decay is ap- 
proximately exp [— (i/T 2 ) 2 ], while the echo can be fitted to exp [— (2t/T 2 ) 2 - 3 ]). 
The problem lies in the fact that the t — > oo assumption averages out finite 
frequency fluctuations; note that T 2 * differs from Ti only via static noise [5*(0) 
in Eq. ([7])]. A large number of pulse spin resonance experiments are sensitive 
to finite frequencies only [the most notable example is the spin echo, which 
is able to remove S z (0) completely]. Therefore one must develop a theory 
for coherent evolution that includes finite frequency fluctuations. Analytical 
results can be derived if we have i) q = for q — x, y (the pure dephasing 
limit) and f] z is distributed according to Gaussian statistics. For many realis- 
tic problems the pure dephasing limit turns out to be a good approximation 
to describe phase relaxation. A common situation is that the neglected com- 
ponents S x ^y{uj) are much smaller than S z (u>) for positive frequencies much 
smaller than the qubit energy splitting. For example, in the case of a localized 
electron spin in a semiconductor, S XiV (u>) cx to 5 due to the combination of the 
spin-orbit and electron phonon interactions [10]. As a result, S Xty (w) <§C S z (ui) 
as lo — > 0. The Gaussian approximation is described below. 

2.2 Finite frequency phase fluctuations and coherence decay in 
the Gaussian approximation 

In many cases of interest, the environmental variable rf z is a sum over several 
dynamical degrees of freedom, and measurement outcomes for the operator 
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f\ z may assume a continuum of values between — oo and +00 In those situ- 
ations we often can resort to the Central limit theorem which states that the 
statistics for outcomes 77' follows a Gaussian distribution, 

^W=^ = ^«p(-^), (8) 

with a stationary (time independent) variance given by A 2 = {fjl(t)). Here 
(A) = Tr{/5eA} is a thermal average taken over all bath degrees of freedom 
(pb oc c- n B/k B T j s ^ ne canon i ca i density matrix for the bath). We also 
assume (fj z (t)) — 0, since any constant drift in the noise can be incorporated 
in the effective B field. 

Our problem is greatly simplified if we relate the operator f] z to a Gaussian 
stochastic process r/(t) in the following way. For each time t, rf{t) corresponds 
to a classical random variable, that can be interpreted as the outcome of 
measurements performed by the qubit on the environment. This allows us 
take averages over the bath states using Eq. {BJ. Note that the statement 
"Gaussian" noise refers specifically to the distribution of noise amplitudes, 
that is not necessarily related to the spectrum of fluctuations (see below). 

Our simplified effective Hamiltonian leads to the following evolution op- 
erator [recall that we set fj x ^ y = in Eq. (JTJ) ] 

H„(t,0) - e-*/o dt ' u W = e -i^[Bt+J^'(t')dt'] = e -4^[s i+ x„(t)] ; (9) 
where we define 

x v (t) = I v'(t')dt'. (10) 
Jo 

Here the subscript r\ emphasizes that this operator is a functional of the tra- 
jectory rj'{t'). The effect of the distribution of trajectories rf can be described 
by assuming the qubit evolves according to the density matri>H 

p{t) = 53^(t,0)poZ<j(t,0), ( n ) 
n 

where p v denotes the appropriate weight probability for each environmental 
trajectory and po is the t = density matrix for the qubit. The coherence 
envelope at time t averaged over all possible noise trajectories is then 

{(a + (t)))=Tr{a + p(t)} 

= ^2 Pr ,Tr{ul(t,0)a+Ur l (t,0)po} 
n 



1 An important exception is the observation of individual random telegraph noise 
fluctuators in nanostructures (in this case i) = ±77' assumes only two discrete val- 
ues). This results in important non-Gaussian features in qubit evolution. 

2 This assumption is equivalent to the Kraus representation in the theory of open 
quantum systems. 
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n 

= (e iX *W)e iBt Tr{a + po}, (12) 

where we used the identity e larT * a + e~ la<7z = e 2la a + . Here the double average 
((■)) denotes a quantum mechanical average over the qubit basis plus an en- 
semble average over the noise trajectories T)'(t). We can evaluate the coherence 
amplitude explicitly by noting that the random variable X(t) — J* r]'(t')dt' is 
also described by a Gaussian distribution, but with a time dependent variance 
given by a t = (X 2 (t)). Therefore we have 

/OO -| _ x 2 

-^^e ^e lX =e-^, (13) 

with at given by 

dt x / dt 2 (v'(tiW(h)) 
'o Jo 

= 2 f dtx [ 1 dt 2 (ri'(t 1 )r]'(t2)) 
Jo Jo 

rt pt-t' /2 

= 2 dt' dT( V '(T + t'/2)rj'(T-t'/2)) 

JO Jt'/2 

= 2 f (t-t')S(t')dt'. (14) 
Jo 

Here we introduced the time dependent correlation function S(t') = (rj'(T + 
t'/2)rj'(T - t'/2)) = (r)(t')r)(Q)), that is independent of a time translation T 
by virtue of the stationarity assumption. Based on the discussion above, it is 
natural to include quantum noise effects by making the identification 

S(t) = ( V '(tW(0)) -> (%(t)7fe(0)). (15) 

It is straightforward to generalize Eqs. (fT3|) and (fT4|) for echo decay In- 
stead of free evolution (denoted free induction decay in magnetic resonance) , 
consider the Hahn echo given by the sequence it/2 — t — n — t — echo. Here 
the notation "7r/2" denotes a perfect, instantaneous 90° spin rotation around 
the x-axis (described by the operator e~i' K(Jl ). The notation "— r— " means 
the spin is allowed to evolve freely for a time interval r. 'V denotes a 180° 
spin rotation around the x-axis, also referred as a "7r-pulse" (this is described 
by the operator e - ^ "* — — io x ). The initial ir/2 pulse prepares the qubit in 
the state po = \y+)(y + |, after which it is allowed to evolve freely for time 
r, when the 7r-pulse is applied. After this pulse the qubit is allowed to evolve 
for a time interval r again, after which the coherence echo is recorded. Hence 
the evolution operator is given by 

W H ahn(2r) =U(2T,T)(-ia x )U(T,0). (16) 
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The same procedure leading to Eq. (Q3J) is now repeated in order to calculate 
the magnitude of the Hahn echo envelope at t = 2r. The quantum average is 
given by 

(o+(2r)) = Tv{U^T,0){ia x )UH2T,T)a + U{2T,T)(-ia x )U(T,0)p } 
= Tr |wt( Tj ) a jC dt ' [B+r '' {t ' )] a + a x U(r, 0)p } 

= Tr^ST dt ' [B+v ' {t ' )] U\T 1 Q)<jM{T,Q) Po 

= e i^dt'[B + r,\t')] e -ip o dt'[B +v '{t')]^ {a _ pQ} _ 

Therefore the double average can be conveniently written as 

«o-+(2t)» =e- iB lT s VW (e-ifr'VMW^Tria-po}, (18) 

with the introduction of an auxiliary echo function s(t). For Hahn echo s(t) = 
1 if < t < t and s(t) = — 1 if t > t. Note that the first term in Eq. (fJU) is 
exactly equal to one, because the Hahn echo is able to completely refocus a 
constant magnetic field. It is convenient to introduce the noise spectrum in 
Eq. (fTgj) via S(t) — J e~ tuJt S(uj)duj in order to get the following expression 
for the coherence envelope 



|<K(t)»|=exp 



du> S , (w)J r (t, lo) 



(19) 



Here we define a filter function that depends on the echo sequence s(t'), 

F(t,w)= I dt's(t') [ dt"s(t")cos[u(t' -t")]. (20) 
./o Jo 

For free induction decay [s(t) = 1] we have 

_ , . lsin 2 (^/2) , \ 

F FID (t,u) = - \ L > , (21) 
2 (w/2) 

while for the Hahn echo 

1 sin 4 (lot/2) 

2 (uj/Af 

Note that ^ r Hahn(2T, 0) = 0. The Hahn echo filters out terms proportional 
to S(Q) in qubit evolution, this is equivalent to the well known removal of 
inhomogeneous broadening by the echo. Any spin resonance sequence con- 
taining instantaneous 7r/2 or 7r-pulses (not necessarily equally spaced) can be 
mapped into an appropriate echo function s(t). An important example is the 
class of Carr-Purcell sequences that can be used to enhance coherence. 



JW2t,o,) = ± OU ; T - (22) 
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General results for the short time behavior 

We may derive interesting results when the time-dependent correlation func- 
tion S(t) is analytic at t = 0. In this case we may expand around t = to 
obtain 

/oo 
-OC 

with the 2n-th moment of the noise spectrum defined as 

/oo 
dujS(uj)w 2n . (24) 

Hence if S(t) is analytic at t = 0, we must have M<i n < oo for all n, i.e. the 
noise spectrum has a well defined high frequency cut-off. It is important to 
keep in mind that the assumption of analyticity at t = is actually quite 
restrictive. Physically, only Mo < oo is required, so that 5(0) < oo (this is 
the noise power or mean square deviation for 77'). Important examples where 
S(t) is not analytic at t = include the Gauss-Markov model described below 
[See Eq. 

If the t = expansion exists we may immediately obtain the short time 
behaviors for the free induction decay and Hahn echo: 

«° + (t)»Fm = e -S duS ^ ! ^r^ 1 « e-!^ot 2 +AM 2 ^ (25) 
« CT+ (2r))) Hah „ = e-^^^W 2 „ e-iM^M^. (2g) 

The short time behavior described by Eqs. (|25|) . (|26p is universal for noise 
spectra possessing a high frequency cut-off. Note the striking difference in 
time dependence: For free induction decay the coherence behaves as ~ e _t , 
while for a Hahn echo we have ~ e~ T . This happens because the Hahn echo 
is independent of the mean square deviation Mq = 5(0). 



Example: The Gauss-Markov model 

The simplest model of Brownian motion assumes a phenomenological corre- 
lation function that decays exponentially in time, 

S(t) = A 2 exp(-\t\/r c ), (27) 

where r c is a correlation time that describes the "memory" of the environ- 
mental noised This model is useful e.g. in liquid state NMR in order to 

3 Many authors use the terminology "Markovian dynamics" to denote evolution 
without memory, i.e., the limit r c — > in Eq. I|27p . This limit can be taken by 
setting A —+ 00 with F = A 2 t c held finite. In that case we have S(t) — > 2r8(t) 
resulting in a "white noise" spectrum and (cr+) oc e~ rt . 
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calculate the line-widths of a molecule diffusing across an inhomogeneous 
magnetic field. In that case, A becomes the typical field inhomogeneity, while 
the "speed" for diffusion is of the order of A/r c . The resulting environmental 
noise spectrum [Fourier transform of Eq. (|27|l ] is a Lorentzian, given by 

A 2 t 1 

S(")= ° ■ (28) 

(UT C ) + 1 



We start by discussing free induction decay. Using Eqs. (|18|) and (|27|) with 
s(t') = 1, we get 

I«^+(*)»Ifid = exp {-Z\ 2 r c 2 [i/r c + (e"*^ - l)] }. (29) 
For t > r c , Eq. leads to 

|((a + (i)))| FID «e-^*. (30) 

In this regime, the correlation function Eq. (|27p can be approximated by a 
delta function, and the decay is a simple exponential signaling that a Master 
equation approach is appropriate [Eq. ([2])], The coherence time is given by 
T 2 * = 1/(A 2 t c ). Interestingly, as r c — > with Zi finite, T 2 * — > oo. This phe- 
nomenon is known as motional narrowing, inspired by the motion of molecules 
in a field gradient. The faster the molecule is diffusing, the narrower is its 
resonance line. Now we look at the low frequency noise limit, t -C r c . This 
leads to 

|((a + (t)»| FID «e-^ 2t2 ^e-^*) 2 . (31) 



In contrast to Eq. (|30|) . the decay differs from a simple exponential and is 
independent of r c . This result is equivalent to an average over an ensemble 
of qubits at any specific time t (in other words, the line- width A/y/2 and de- 
phasing time T 2 = y2/A are a consequence of inhomogeneous broadening). 
Therefore the coherence decay is completely independent of the environmen- 
tal kinetics. As we shall see below, this decay is to a large extent reversible 
by the Hahn echo. 

The Hahn echo decay is calculated from Eqs. (flU)) . (f2"2")) , and (|2"5|) leading 

to 

«o+(2T)»Hahn = exp {-A 2 t 2 c [2t/t c - 3 + 4c~^ - e" 2 ^] }. (32) 

For r 3> t c we again have motional narrowing, ((cr_|_(2r)))Hahn ~ e ~ A t c 2t j a 
result identical to FID [Eq. (f3"0"|) ] if we set t = 2r. This occurs because the 
noise trajectories are completely uncorrelated before and after the 7r-pulse, 
that plays no role in this limit. For t<t c we get 

((a + (2r))) H ahn w e~^^) 2 (^) 3 « e-(M)\ (33) 
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In drastic contrast to free induction decay, Eq. (|33p depends crucially on 
the kinetic variable r c . The time scale T 2 for 1/e decay of Hahn echo^ is 
considerably longer than T 2 when A^$> t c . 



A train of Hahn echoes: The Carr-Purcell sequence and coherence 
control 

Consider the sequence 7r/2 — [t — tt — t ~ e c ho] repoat . It consists of the appli- 
cation of a 7r-pulse every odd multiple of r, with the observation of an echo 
at even multiples of t, i.e. at t — 2nr for n integer^ In the limit r <C t c the 
n-th echo envelope can be approximated by a product of n Hahn echoes, 

«a + (2nr))) CP « «a + (2r))}£ ahn « e^i^T = e~^, (34) 

with T| ff = T 2 [T 2 /(2r)] 2 . As r is decreased below T 2 the effective coherence 
time T| increases proportional to 1/r 2 . Therefore a train of Hahn echoes 
can be used to control decoherence. Rewriting Eq. (|34[) with t = 2nr we get 
T| ff = (2n) 2 / 3 T 2 , showing that the scaling of the enhanced coherence time 
with the number of 7r-pulses is sub-linear. The train of 7r-pulses spaced by 
r -C t c effectively averages out the noise, because within r much shorter than 
t c the noise appears to be time- independent. 



Loss of visibility due to high frequency noise 



In order to understand the role of high frequency noise, consider the model 
Lorentzian noise spectrum peaked at frequency Q with a broadening given 
by 1/T d , 

A 2 r d 1 



Sl(uj) — 

« (co-fifri + l 

Using Eqs. p^]) and (|2~lj) and assuming J? ^> 1/r^ we get 

2 



exp 



1 - e~ t/Td cos Qt 



(35) 



(36) 



Therefore high frequency noise leads to loss of visibility for the coherence 
oscillations. The loss of visibility is initially oscillatory, but decays exponen- 
tially to a fixed contrast for t ^> t^. For comparison consider the Gaussian 
model, 



4 In the electron spin resonance literature the 1/e decay time of a Hahn echo is 
often denoted Tm- Here follow the spintronics terminology and use T2 for the 1/e 
decay of Hahn echo, and for 1/e decay of FID in the low frequency regime. 

5 One can make the Carr-Purcell sequence robust against pulse errors by alter- 
nating the phase of the 7r-pulses, see e.g. the Carr-Purcell-Meiboom-Gill sequence 
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2a 2 



(37) 




) 



2 



l«^+W»l G «exp 



( 



1 — e 2 



cos fit 



) 



(38) 



Note that the difference between the Gaussian and the Lorentzian models lies 
in the time dependence of the approach to a fixed contrast. Deviations from 
Lorentzian behavior may be assigned to non-exponential decays of the co- 
herence envelope. Although Eqs. (|3l)|) . (f3"8"]) were calculated for free induction 
decay, it is also a good approximation for Hahn echoes in the limit Q ^> 1/r. 

2.3 Single spin measurement versus ensemble experiments: 
Different coherence times? 

Recently single shot detection of the spin of a single electron in a GaAs quan- 
tum dot was demonstrated [4] , and the Hahn echo decay of the singlet-triplet 
transition in a double quantum dot was measured [7]. Also, spin resonance 
of a single spin center in the Si/Si02 interface was detected through time 
averaged current fluctuations [5]. These state of the art experiments should 
be contrasted with the traditional spin resonance measurements where the 
microwave excitation of a sample containing a large number of localized spins 
is probed. Naturally the following question arises: Are the coherence times ex- 
tracted from an ensemble measurement any different from the ones obtained 
in single spin experiments? 

The answer to this question is related to the ergodicity of the environment 
producing noise, i.e., whether time averages are equal to ensemble averages. 
Even when single shot read-out of a quantum degree of freedom is possible, 
one must repeat each measurement several times, in order to obtain good 
average values for the observables. For example, measurements of the state 
of a single spin yields two possible outcomes and one must time-average an 
ensemble of identical qubit evolutions in order to obtain an average value 
that reflects the correct outcome probabilities (In Ref. [I] , each average value 
resulted from ~ 600 read-out traces). The presence of phase fluctuation with 
correlation time r c smaller than the typical averaging time implies spin pre- 
cession with distinct frequencies for read-out traces separated in time. This 
may lead to strong free induction decay (low T 2 *) m a single shot read-out 
measurement, see Fig. [2] This was indeed observed in the double-dot exper- 
iments of Petta et al [7]. 

The free induction decay time T£ in ensemble experiments may be quite 
different from single-spin experiments. This is because the spatial separation 
of spins adds several new contributions to zero frequency noise. These include 
spatially inhomogeneous magnetic fields, g-factors, and strains. However, all 
these contributions are removed by a Hahn echo sequence, that filters out 



Loise 13 



(a) Ensemble of spatially separated qubits 




(b) Time average of a single qubit 




Fig. 2. (a) A Traditional spin resonance experiment probes the coherent evolution 
of an ensemble of spatially separated spins. The spins can be separated into different 
"packets" with similar resonance frequencies, each packet with linewidth 2/T2. A 
free induction decay measurement is sensitive to the broadened linewidth 2/T2 . A 
spin echo is needed in order to reveal the intrinsic linewidth 2/T2. (b) A similar 
situation applies to single spin experiments subject to low frequency noise, because 
each time ensemble may have a different resonance frequency 



environmental fluctuations with frequencies smaller than 1/r. Therefore we 
expect the spin echo decay time T2 for a spatial ensemble and a single spin 
to be similar, provided the mechanisms for finite frequency noise do not 
vary wildly for spatially separated spins, and the environment affecting each 
individual spin is ergodic within the average time scale (similar reasoning 
applies for T± , which is determined by high frequency noise uj = 7S) . 

The Gaussian theory of decoherence described here is appropriate for 
ergodic environments. It would be very interesting to explore model systems 
both experimentally and theoretically in order to search for detectable non- 
ergodic effects in coherent evolution. 



3 Electron spin evolution due to nuclear spins: Isotropic 
and anisotropic hyperfine interactions, inter-nuclear 
couplings and the secular approximation 

A localized electron spin coupled to a lattice of interacting nuclear spins 
provides a suitable model system for the microscopic description of environ- 
mental fluctuations affecting coherent evolution. Here we describe a model 
appropriate for localized electron spins in semiconductors, and discuss some 
approximations that can be made in a moderate magnetic field (typically 
larger than the inhomogeneous broadening linewidth, B > 0.1 T). 
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3.1 The electron-nuclear spin Hamiltonian 

The full Hamiltonian for a single electron interacting with N nuclear spins is 
given by [23] 

H = H e Z + HnZ + W en + H en + ~Hnn- (39) 

Here the Zeeman energies for electron and nuclear spins are respectively 

HeZ = \leBo z (40) 

H n z = -7nBj2 I i^ ( 41 ) 

i 

where ar = (a x , o y , a z ) is the Pauli matrix vector representing the electron 
spin, and Ii is the nuclear spin operator for a nucleus located at position .R, 
with respect to the center of the electron wave function. For B = 1 T we 
have 7e-B ~ 10 11 Hz while j n B ~ 10 s Hz. The e — n coupling takes place due 
to isotropic and anisotropic hyperfine interactions. The isotropic hypcrfinc 
interaction is given by 

W e „ = ^Afli-^, (42) 

i 

with contact hyperfine interaction 

4 so = Y^nhMRi)\ 2 , (43) 

where j e o = e/{m e c) = 1.76 x 10 7 (sG) _1 is the gyromagnetic ratio for a 
free electron and W(r) is the electron's wave function. Typical values of Af° 
varies from Af ° ~ 10 7 Hz for Hi = (at the center of a donor impurity wave 
function) to Af° ~ for i?, much larger than the impurity Bohr radius. The 
anisotropic hyperfine interaction reads 

W'en = ~I>-^.<7 (44) 

i 

with an anisotropic hyperfine tensor given by 

{A' i)lm = le0ln H J *rMr ^ (^) - ^) , (45) 



where = e 2 /(m e c 2 ) is the classical electron radius. Note that Eqs. ((421 
and (HU) are first order perturbative corrections in the electron coordinate r 
Finally, the nuclear-nuclear dipolar coupling reads 



Ii ■ Ij 3 (Ii ■ R i3 ) (Ij ■ Rij) 



R 3 R 5 - 



(46) 

where Rij = Ri — Rj is the distance between two nuclei. The typical energy 
scale for Eq. (|16")) is a few KHz for nearest-neighbors in a crystal. 
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The full Hamiltonian Eq. (|39p is a formidable many-body problem. It 
is particularly hard to study because of the lack of symmetry. In order to 
study theoretically the quantum dynamics of an electron subject to a large 
number of nuclear spins we need to truncate Eq. (|39[) . Here we discuss some 
simplifications appropriate for JS > 0.1 T, a condition typically satisfied in 
several experiments. The first approximation arises when we note that the 
electron Zeeman energy is typically 10 3 times larger than the nuclear Zeeman 
energy. For B > 100 G the former is much larger than Af°, therefore the 
electron spin cannot be "flipped" by the action of the hyperfine interaction. 
In other words, "real" e — n flip-flop transitions get inhibited at these fields 
(however virtual transitions induced by second order processes such as Hl n 
do produce visible effects as discussed in section I3.3[) . This consideration 
allows us to approximate the isotropic hyperfine interaction to a diagonal 
form (secular approximation), 

«en«^£>ri i2 . (47) 

i 

The anisotropic hyperfine interaction contains a similar diagonal contribution 
in addition to pseudo-secular terms of the form a z Ii±. These terms lead 
to important echo modulations of the order of ~ (A^/^fnB) 2 ~ 0.1 — 1 for 
moderate magnetic fields (B ~ 0.1 — 1 T). To derive these terms, assume b as 
the direction of the magnetic field and substitute <x — > a z b in Eq. (|4"4")) . The 
result is 

H' en « \o z (A' q I tz + A' l± I l+ + A'*Ji.) , (48) 
where 

A' q = z ■ A\ ■ b, (49) 



Al I 



\ [x ■ A\ ■ b + ly ■ A' r 1] . (50) 

For some lattice sites (closer to the center of the donor impurity) A' i± is a 
reasonable fraction of the nuclear spin Zeeman energy (even for B ~ 1 Tesla), 
and as a consequence the precession axis of these nuclear spins will depend 
on the state of the electron, producing strong modulations in the nuclear spin 
echo signal For the electron spin Eq. (jig]) produces small modulations 
observed at the shortest time scales in the echo decay envelope [19]. Finally, 
we can neglect contributions to Eq. (|46|) that do not conserve nuclear spin 
Zeeman energy, 

n nn w V hi (ii+ij- + - Mizijz) , (5i) 



i<j 



with 



1 ^1-3C0S 2 ^; 



ba = ~4< n ^3 (52) 
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Here 6^ is the angle formed by the applied B field and the vector i£y linking 
the two nuclear spins This leads to an important orientation dependence 
of coherence times. 

3.2 Electron-nuclear spin evolution in the secular approximation 

In the secular approximation [Eqs. ([47]) - ([51]) ] the electron-nuclear spin Hamil- 
tonian is block-diagonal, 

H = H + |t)(t|+H_||)(4|, (53) 

where and |-U-)(-U-| are projectors in the electron spin up and down 

subspaces respectively. Here TL± contains only nuclear spin operators and is 
given by 

1 

i 

i 

where Ai = Af° + AL and A\ = A' i± . Accordingly, the evolution operator 

becomes U{t) = U+(t) +U-(i) |JJ-)(JJ-|, with U±{t) = e" jm ±. We can 

write an explicit expression for the coherences as a function of the evolution 
operators U± (t) provided the initial density matrix can be written in product 
form, po = poe <8> pa n - The free induction decay is given by 

(cr+(t)) F iD = Tr n Tr e {U f a + Upa c (g) p 0n } 

= Tr n I^IW-ponj Tr c {<J + p 0c } , (55) 

where we used the fact that U f a + U = u\.U- |lr)(JJ-|, and a+ =|lr)(# For the 
Hahn echo, we use U(t)cj x U(t) as our evolution operator, to get 

(a + (2r)) H ah„ - Tr n \^_u\UM+p Qn } Tr c {<T- POc } . (56) 

Eqs. (|55p and (|56p are exact in the secular approximation, and make ex- 
plicit the dependence of the electron's coherence envelope in the nuclear spin 
Hamiltonian evolution. 

Inhomogeneous broadening due to the isotropic hyperfine 
interaction 

The diagonal model 

i i 

is easily solved exactly for nuclear spins initially in a product state. Assume 
the electron spin is pointing in the +y direction, and the nuclear spin states 
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are distributed randomly, each nuclei with equal probability of pointing up 
or down. The free induction decay amplitude becomes 



II 5^ 



i e iy*LH ]"f 1 \ r -.\j ,- v,/ 



*vy.m \- J_„ §*E 1 « 

fi=±b- 



~ l e iTeBt e -s* 2 E^4 ; ( 58 ) 

where in the last line we assumed AT — > oo with each individual Aj — > 
so that the hyperfine field can be approximated by a continuous Gaussian 
distribution. The free induction decay rate or inhomogeneously broadened 



linewidth is given by ^ ~ j4 rms = y Ej A*. This fast decay rate should be 

compared to the Hahn echo: From Eq. (|56p we see that (cr+(2r))Haim = — i/2, 
the Hahn echo never decays. In fact, from Eq. (|56[) we can easily prove that 
the class of Hamiltonians of the form Eq. ([55)1 satisfying [7i + ,7Y_] = have 
time independent Hahn echoes given by Tr{cr_p 0e } [2"5] , 



3.3 Beyond the secular approximation: Nuclear-nuclear 
interactions mediated by the electron spin hyperfine interaction 

In the sections above we showed that the secular approximation allows us 
to decouple electron spin dynamics from nuclear spin dynamics completely. 
This approximation clearly does not hold at low magnetic fields, and the 
problem becomes considerably more complicated. The study of electron spin 
evolution subject to the full isotropic hyperfine interaction has attracted a 
great deal of attention lately [551 [57J [551 [55] , particularly because of a series of 
free induction decay experiments probing electron spin dynamics in quantum 
dots in the low magnetic field regime [29) . In the author's opinion the most 
successful theoretical approach so far in the description of these experiments 
is to treat the collective nuclear spin field classically by taking averages over 
its direction and magnitude [55]. Here we shall not discuss the interesting 
effects occurring at low fields. Instead, we will focus on the following question: 
What is the threshold field B t h for the secular approximation to hold? For 
intermediate B > B t h (not satisfying B ^> Bth), how can the non-secular 
terms be incorporated in a block diagonal Hamiltonian of the form Eq. (|53l) ? 
In order to answer these questions, let's consider the Hamiltonian 

H = H Q +V, (59) 



= ^(Te + ln)Ba z + -a z Ajl jz , 

j 



(60) 
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V 



j 



A 3 {a+Ii- 



<J-I ]+ ) . 



(61) 



where Ho and V denote the secular and non-secular contributions respec- 
tively. Here we remove the nuclear Zeeman energy by transforming to the 
rotating frame precessing at j n B. From Eq. (fSTj) we may be tempted to as- 
sume that flip-flop processes involving an electron and a nuclear spin (e.g., 
ft! - MJ-t) are forbidden by energy conservation at fields j e B 3> A4. However, 
the situation is much more complex because higher order "virtual" processes 
such as V 2 = AiAjIi+Ij_a+cr- + ■ ■ ■ preserve the electron spin polarization 
and hence may have a small energy cost (of the order of Aj — Aa for Ai ~ Aj ) . 
As we show below, these processes actually lead to a long-range effective cou- 
pling between nuclear spins, similar to the RKKY interaction between nuclear 
spins in a metal. We will show this using the original self-consistent approach 
of Shenvi et al. 

Let be a "+" eigenstate of the Hamiltonian Eq. (l59|) . i.e. has 
primarily electron spin-up character. Without loss of generality, \ip + ) can be 
written as 



(62) 



Because the perturbation V flips the polarization of the electron, the action 
of H on the electron spin-up and electron spin-down subspaces yields the two 
simultaneous equations, 



Ho 
Ho 



Eq. lHHD can be solved for 
Eq. (631) yields 



= £4 



(63) 
(64) 



and the resulting expression inserted into 



Ho 



V- 



1 



E+~Ho 



-V 



E 4 



(65) 



In the presence of an energy gap between the spin-up and spin-down states 
(this is certainly true at high magnetic fields satisfying B > Aj Me), the 
operator 1/(E+ — Ho) is always well-defined [30]. Because the left-hand side 
of Eq. (|65|) depends on E+, it is not a true Schrodinger equation; to obtain 
E+ exactly, Eq. (|65| must be solved self-consistently. However, if we use 
E + « (7 e + j n )B/2, then we can obtain an effective Hamiltonian from Eq. 
(|65[) . The effective Hamiltonian in the electron spin-up subspace is 



Hi 



Ho 
1 V - 



v e + ff 



3,k 



1 



( 7e + ln )B 



2 Ajljz 



h.4 



(66) 
(67) 
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We obtain a similar, but not identical, effective Hamiltonian for the spin-down 
subspace (note the transposition of the 7_ and I + operators), 

«i = Ho + V,j (68) 

Eqs. ([6"T]) and ([6"9"]> show that the overall coupling between nuclei does indeed 
decrease at high fields, because the operator l/(E — Tto) scales approximately 
as 1/B. However, the energy cost for flip-flopping two nuclei j and k is pro- 
portional to Aj — Ak- Thus, if Aj and Ak are close in value, the nuclei can 
flip-flop even at high fields. Eqs. ([67]) . ([69]) were later derived using an alter- 
native canonical transformation approach [3"T] . 

We may expand Eqs. ([67)) and ([69]) in powers of J^j AjIj Z / (7e-B), so that 
for the unpolarized case we have approximately 

1 v-^ AjA k , . 

v— r ^2^3B li+h - <70) 

This effective Hamiltonian is of the secular type Eq. ([53"]) , and satisfies the 
symmetry condition [H+,Ti.-] = 0. Therefore a Hahn echo is able to refocus 
this interaction completely. 

Shenvi et al performed exact numerical calculations of electron spin echo 
dynamics in clusters of N — 13 nuclear spins [2S]. The Hahn echo envelope 
was found to decay fast to a loss of contrast given by 

K° + (2t))1*1- - EjA L 2 - (71) 



2 



[{ le+ln )B} 



This shows that the threshold field for neglecting the non-secular isotropic hy- 
perfine interaction in the Hahn echo is given by the inhomogeneously broad- 



ened linewidth, B t h — A* he ^ 10 — 100 G (for a donor impurity 

in silicon, i? t h is w 1 G for natural samples and ss 10 G for 29 Si enriched 
samples). 

Recently, Yao et al [31] and Deng et al. [28] showed that the electron- 
mediated inter-nuclear coupling may be observed as a magnetic field de- 
pendence of the free induction decay time in small quantum dots. There is 
currently an interesting debate on the correct form of the time dependence 
for FID decay. Yao et al. derived the FID decay from Eq. ([70]) and obtained 
((7-1- (i)) ~ e~* , while Deng et al. carried out a full many-body calculation to 
argue that FID scales as a power law according to (a + {t)) ~ 1/t 2 . 
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4 Microscopic calculation of the nuclear spin noise 
spectrum and electron spin decoherence 

In this section we discuss low frequency noise due to interacting nuclear spins. 
For simplicity, we assume only isotropic hyperfine interaction. The inclusion 
of anisotropic hyperfine interaction is considerably more complicated, but can 
be seen to lead to high frequency noise and echo modulations (at frequencies 
close to 7 n i? ~ 10 s s^ 1 per Tesla). In the next section we provide explicit 
numerical calculations for the case of a phosphorus impurity in silicon and 
compare our results to experiments. 



4.1 Nuclear spin noise 

Using the approximations Eqs. (|47|) and ([51]) we can write the electron-nuclear 
Hamiltonian in a form similar to Eq. {!]). In the electron spin Hilbert space, 
we assume an effective time-dependent Hamiltonian of the form 

Ks = \leBa z + l -a z J2 Mhz (*), (72) 

i 

with Ai = Af°. Nuclear spin noise is in turn determined by the effective 
Hamiltonian 

Tiij = -f n B (I lz + I jz ) + b i:j (Ii+Ij- + h-Ij+ - ^hzljz) 

+i {Ail iz + Ajl jz ) , (74) 

where we decoupled the electron spin from the nuclear spins by assuming the 
nuclear spin wave function evolves in the electron spin up subspace (a z — > 1). 
An equally valid choice is to assume a z — ► — 1. It turns out that this choice 
does not matter within the pair approximation described below. We will 
check this by noting that the final answer is unchanged under the operation 
Ai -Ai for all iE 

For now we assume the nuclear spins are unpolarized (T = oo) so that 
Aili z ) = 0. This approximation will be relaxed below. The time depen- 
dent correlation function for nuclear spins is given by 



As we shall see below, this approximation leads to identical results as the lowest 
order cluster expansion developed in [21] . However, interesting interference effects 
arise when this approximation is not valid. The cluster expansion method beyond 
lowest order [22] takes account of the full electron-nuclear evolution, therefore it 
can be used to study these effects. 



t„b 




TT 



|+} = aU)+/?lT 
\-) = -fi\1l)+a\lT) 



Fig. 3. Energy levels for two nuclear spins coupled through the dipolar interaction. 
The flip-flop mechanism corresponds to transitions between the states |+) and | — ), 
which are admixtures of and |J.| states. The anisotropic hyperfine interaction 
couples states differing by ~ 7n-B in energy. 



S(t) = (y^Aili^^AjIjM) 

= ^4 2 </, z (^(0))+ Yl AA 3 (I lz (t)I 3Z (0)). (75) 

We now invoke a "pair approximation" by assuming 

(/te(t)/i»(0)) « X>i*(*)J«(0))y, (76) 

(I iz (t)I jz (0)) « (Q)>0, (77) 

where (-)y denotes a thermal average restricted to the ij Hilbert space. The 
operator Ii Z (t) is in the Heisenberg representation defined by the two particle 
Hamiltonian Eq. (74]). Plugging Eqs. (76]) and (77]) in Eq. ([75]) and reordering 
terms we get 



with 



(78) 



(79) 



The same derivation could be given for finite temperature, when the thermal 
average of the hyperfine field is non-zero. The result is identical to Eq. (|78p 
except for the substitution fjij — > Sfjij = fjij — {fjij)- Using the definition 
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of the noise spectrum [Eq. ©] and expanding the correlator in the energy 
eigenstates of the pair Hamiltonian Eq. (TT4"1) we get 

S l3 (to) = J2p a \(a\Sm J \(3)\ 2 5 (lu - E Pa ) , (80) 

a, 13 

with Epos = Eft — E a the energy difference between the energy eigenstates 
| a), andp Q the (thermal) occupation of state a. Therefore the noise spec- 
trum is a sum over all possible transition frequencies induced by the operator 
For nuclear spin 1/2 the ij Hamiltonian has the following eigenenergies 
and eigenstates (See Fig. [3]) 

= ~f n B - by + <Hj, (81) 



E + = b l3 + Jbl+A%, (82) 



E-=b ij - s ]bl+A%, (83) 
E u =--y n B-b ij -a ij , (84) 



|+> = cos- ItD + sin- lit), (85) 
|-> = -sh4 |U)+co4 UT), (86) 



with 



2 2 

di^^iAi+Aj), (87) 

Aij =\{Ai- A,) , (88) 

cos 6> = — p =l (89) 



b 2 +A 2 . 

K 

b 2 - + A 2 - 



(90) 



Using Eqs. {55} , (|86p the transition matrix element is easily found to be 

6 

H»tol+) = - 4 Aj sin - cos - = -24 y sin0. (91) 
The transition frequency is simply the difference between Eqs. (|82|) and (|83| . 



£+_ = 2^bf + ^.. (92) 
The resulting noise spectrum is therefore 

Sij(w) = (A™ s ) 2 S(cj)+4 + £+_) + - £+_)] , (93) 

with a static contribution given by 



(-4 



rms\ 
ij ) 
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At 



(ptt + Pu) a « + (JP++P-) 



b 2 +A 2 - 



(PTT -Pii) a H + (p+-p-) 



A 2 - 



b 2 +A 2 - 



(94) 



The free induction decay due to this noise spectrum can be easily calcu- 
lated using Eq. (fT9|) and the filter function Eq. (|2"Tj) . 



«o + (t)» = exp 



x exp 



2 ^ 



)V 



-2E(p-t 



6 2 - 



Z\ 2 . 



< 2 sinc 2 



(v^ 



,(95) 



where sinc(x) = sin(x)/x. As expected, the FID decay is usually dominated 
by the zero-frequency noise amplitue A™ [Eq. (|94|) ]. The dipolar- induced 
decay may be visible provided A™ s is much smaller than the finite frequency 
noise amplitudes. For example, this is the case if the nuclear spins are polar- 
ized, since we have A™ s = exactly when = 1 in Eq. (I94|) . 

The Hahn echo decay envelope is derived after integration with the filter 
function Eq. f22l. 



«ct + (2t))) =exp 



-8Y,(P++P-)b^T i sinc 

i<j 



A 2 -T 

13 



(96) 



At T — > 00 (p+ = = 1/4) Eq. (fM]) is identical to the echo decay obtained 
by two completely different methods, viz. the lowest order cluster expansion 
[Eq. (20) in Ref. [21]] and the quasiparticle excitation model [Eq. (18) in 
Ref. [31J]. As expected, Eq. |96|) is independent of zero-frequency noise, and 
is exactly equal to 1 when either 6 y = or Z\ y = 0, or when the nuclear 
spins arc polarized (p+ = p- =0). 

By expanding the exponent in Eqs. (|95| and (|96| in powers of time, we 
find that only even powers are present. The short time behavior for FID 



is 



while for Hahn echo 



This short time approximation is 



valid for times much shorter than the inverse cut-off frequency of the noise 
spectrum obtained after summing over all pairs i, j. 



4.2 Mean field theory of noise broadening: quasiparticle lifetimes 

We showed that the noise spectrum due to flip-flop transitions in the Hilbert 
space formed by two nuclear spins i,j is a linear combination of delta func- 
tions. We may extend this pair approximation to clusters larger than two, 
and the number of delta functions will grow exponentially with cluster size 
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Fig. 4. An isolated pair of nuclear spins with infinite lifetime will produce a sharp 
peak in the noise spectrum. The role of the many-body interactions with other 
nuclear spins is to broaden this peak and smooth out the noise spectrum for the 
collective nuclear spin excitations. Here we calculate the line broadening for each 
pair flip-flop transition using a procedure similar to van Vleck's method of moments. 



[This can be done by sistematically increasing the size of the Hilbert space 
beyond a single pair i,j in Eqs. (|7rj|) and (|77[) ]. These delta functions can 
be interpreted as well defined nuclear spin excitations with infinite lifetime^ 
With this interpretation in mind, it is natural to expect that the many-body 
interactions of the pair i,j with other nuclear spins will limit the lifetime of 
the flip-flop excitations, i.e., the delta functions will be broadened, see Fig.UJ 
We can add broadening to the delta functions in a mean field fashion by 
using the method of moments, which is applicable at infinite temperature (no 
nuclear spin polarization, i.e. ksT 3> J n B)- In this limit, the noise spectrum 
is written as 

= E ^IH%'I/3)| 2 ^ - fy a ), (97) 

where here a, (3 denote exact many-body eigenstates of the system of N- 
coupled nuclear spins. The n-th moment J uj n S{uj)duj can be calculated ex- 
actly using the invariance of the trace|f| Consider the zeroth-moment, 

Sij(u)du = ^ 22(u\fjij\f3){P\f)ij\a) 



7 In Ref. [3T] Yao et al derive a similar quasiparticle picture via direct calculation 
of the time dependent correlation function for the electron spin. However, the au- 
thors did not calculate the quasiparticle relaxation times. The noise spectrum is a 
natural starting point for developing a theory for quasiparticle energy broadening, 
as we show here. 

8 A similar method was used in the semiclassical theory of spectral diffusion in 
order to calculate the flip-flop rates for pairs of nuclear spins [16| . 



Electron spin as a spectrometer of nuclear spin noise 



25 



Accordingly, the second moment is given by 

1 



A; 



(98) 



lo Sij{uj)doj 



Tr 



(99) 



Note that Eqs. (|98[) and ([99]) are exact at infinite temperature. 

The mean field approximation employed here assumes each delta function 
in the noise spectrum can be approximated by a Gaussian function normalized 
to one0 The noise spectrum becomes 



\<x\mj\ 



at,P 



: exp 



2^ 



and the second moment is 



/°° l 
uj 2 S t:j (uj)duj w ^-|( a |% 



T a(3 



(100) 



(101) 



We now calculate the broadenings <j a p by equating Eq. (jlOip with Eq. (|9"9"| . 
This procedure can be carried out exactly, since the noise spectrum has two 
identical peaks (Fig. 2]) at frequencies ±E + -. The broadening is found to be 



A 2 - 



AA 2 b 2 - 



(blA 2 



(102) 



When A-, <~ b 



■y, the broadening becomes of the same order of magnitude 

b 2 



as van Vleck's second moment for the dipolar interaction (equal to 9 J2k u ik 



[2"5]). For Aij 3> hj we have ct + _ ~ Z\y, and cr + _/i? + _ ~ 1. This type 
of excitation is of high frequency and short lifetime, showing fast decay to 
a small loss of contrast as described by Eqs. (|36| . (|38|) . The broadenings 
describe the diffusion of localized nuclear spin excitations (deviations from 
thermal equilibrium) over length scales greater than the pair distance. 



9 We can verify this assumption by calculating the skewness (fourth moment 
divided by three times the second moment squared). For a perfect Gaussian the 
skewness is exactly one. We carried out this calculation and showed that for Aij <~ 
bij the skewness is very close to one. On the other hand for Aij S> bij the skewness 
becomes large, and a better approximation is a Lorentzian with a cut-off at the 
wings. Nevertheless, by inspecting Eq. (|96[1 we note that nuclear spin pairs with 
Aij ^> bij give a much weaker contribution to echo decay than pairs with Aij ~ bij. 
Therefore this Gaussian fit is precisely valid for most important pairs. As discussed 
in Eqs. (I36|l , (|38p the difference between a Gaussian and a Lorentzian fit lies in the 
time dependence of the decay of coherence modulations; This is ~ e _CT * for a 
Gaussian and ~ e~ t l Td for a Lorentzian. 
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By adding broadenings to the delta functions in Eq. (j9"3"]l we are able to 
plot a smooth noise spectrum, and study the relative contributions of a large 
number of nuclear spins as a function of a continuous frequency. The modified 
Eq. (|93p summed over all nuclear pair contributions reads 



b 2 A 2 



b 2 - 



A 2 . 



P+- 



+ 



i<j 



P-- 



(103) 



For studies of echo decay we may drop the delta function contribution at 
zero frequency. Note that the first part of Eq. (|103[) gives an additional zero 
frequency contribution, that is the limit lu — ► of the broadened spectrum. 



5 Electron spin echo decay of a phosphorus impurity in 
silicon: Comparison with experiment 

In this section we apply our theory to a phosphorus donor impurity in bulk 
silicon. We consider both natural samples (/ = 4.67% 29 Si nuclear spins) 
and isotopically enriched samples (/ = 99.23% 29 Si nuclear spins). We show 
explicit numerical calculations of the nuclear spin noise spectrum resulting 
from dipolar nuclear-nuclear couplings, predict the Hahn echo envelope and 
compare our results with the experimental data of Tyryshkin et al |18j and 
Abe et al [15]. 



5.1 Effective mass model for a phosphorus impurity in silicon 



Here the donor impurity is described within effective mass theory by a Kohn- 
Luttinger wave function [32j . 



1 

^( r ) = ^E^( r K'( r ) 
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with envelope functions Fj describing the effective mass anisotropies. Here 
n = (0.029eV r /£'i) 1 / 2 with Ei being the ionization energy of the impurity 
(Ei = 0.044 eV for the phosphorus impurity, hence n = 0.81 in our case), 
asi = 5.43 A the lattice parameter for Si, a — 25.09 A and b — 14.43 A 
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characteristic lengths for Si hydrogenic impurities [33] . Moreover, we will use 
experimentally measured values for the charge density on each Si lattice site 
|u(Ri)| 2 = 77 « 186 [32]. Hence the isotropic hyperfine interaction is given by 

4*° = ^7e07»^[JMRi)«»(koXi) 

+F 3 (R i ) cos {k Yi) + F 5 (Ri) cos {k Q Zif , (107) 

with the Si conduction band minimum at fco = (0.85)27r/agi, gyromagnetic 
ratios for 29 Si nuclear spins j n = 5.31 x 10 3 (sG)" 1 , and the free electron 
7 e o = 1.76 x 10 7 (s G) _1 . It is instructive to check the experimental validity 
of Eq. (|107p by calculating the inhomogeneous line- width ~ l/(7 e0 T 2 *). A 
simple statistical theory [Eq. ([58]) ] leads to 

(( W / 7e o - Bf) = -L- (4 S °) 2 • (108) 

For natural silicon (nuclear spin fraction / = 0.0467) our calculated root 
mean square line-width is equal to 0.89 G. On the other hand, a simple spin 
resonance scan leads to 2.5 G/2V2 In 2 = 1.06 G [33]. Therefore the simple 
model employed here is able to explain 84% of the experimental hyperfine 
line-width. This is the level of agreement that we should expect when com- 
paring our theory for echo decay with experiment. 



5.2 Explicit calculations of the nuclear spin noise spectrum and 
electron spin echo decay of a phosphorus impurity in silicon 

The nuclear spin noise spectrum is calculated from Eq. (|103[) by excluding 
the 5(u>) contribution. For each pair i, j we calculate the transition frequency 
Eq. (|92|) and broadening Eq. (|102p using the derived microscopic values of 
the hyperfine interaction [Eq. (|107|) ] and the dipolar interaction 

1,1-3 cos 2 0a 

h^-jllh -3 ^ (109) 

For silicon the sites i,j lie in a diamond lattice with parameter a S i — 5.43 A. 
We wrote a computer program that sums over lattice sites Ri within ro of the 
center of the donor. Each site Ri is then summed with all sites Rj within r' 
of Ri (excluding double counting). After numerical tests we concluded that 
the values ro = 200 A and r' Q — 10 A were high enough to guarantee conver- 
gence (increasing tq and r' changes the calculations by a negligible amount). 
Our explicit numerical calculations for the echo decay without broadening 
[Eq. (|Tj6]) ] reproduced the equivalent calculation of Witzel et al. [21] with no 
visible deviation. For ksT 3> ^ n B we may assume that the nuclear spins are 
completely unpolarized (the experimental data was taken at T = 4 K and 
B = 0.3 T [H]). We account for the isotopic fraction / (ratio of sites con- 
taining nuclear spin 1/2) using a simple averaging method. For example, the 
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Fig. 5. Nuclear spin noise affecting the electron spin phase. We show the noise 
spectrum for several magnetic field orientation angles 6 with respect to the (001) 
direction. As 9 is increased from zero, a broad peak develops at a frequency close 
to the dipolar splitting between nearest neighbors in the lattice. In this case the 
noise spectrum is clearly non-monotonous, and cannot be described by a Markovian 
model. The spin echo envelope is a frequency integral of the noise spectrum weighted 
by a filter function. 

pair populations are set as p+ = p_ = / 2 /4, and the broadening <j\_ cx / 
[note J2k b l in Eq. (EE)]. 

Fig. [5] shows the nuclear spin noise spectrum for natural Si at four differ- 
ent magnetic field orientation angles 9 with respect to the crystal direction 
(001). Here 9 = 0° corresponds to B || (001), while 9 = 90° corresponds to 
B || (110). For 9 away from zero the noise spectrum is charaterized by a broad 

peak at which the flip-flop transition frequencies accumulate. The fact 

that the spectrum is non-monotonic implies important non- Markovian be- 
havior for electron spin dynamics (recall that a Markovian noise spectrum is 
defined as a sum of Lorentzians, hence it is always monotonic). Interestingly, 
for 9 close to zero and at low frequencies (cu < 5 x 10 3 s _1 ), the spectrum 
appears to be similar to a Lorentzian peaked at u) = However, one can 
not fit a Lorentzian up to high frequencies because the assymptotic behavior 
deviates significantly from 1/lo 2 . 

The Hahn echo is obtained by integrating the noise spectrum multi- 
plied by the filter function Eq. (ffi?)) up to a frequency cut-off A (we used 
A = 10 6 — 10 7 s _1 , and df ~ 1 — 10 s _1 in our numerical calculations). The 
result is shown in Fig. [S] for two different orientations. We show calculations 
of the echo without broadening [Eq. ([96]) . identical to the result shown in 
Ref. [2Tj] and for the echo with broadening, that is obtained through direct 
integration of the noise spectrum shown in Fig. [5j Note that the two theories 
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t [ms] 

Fig. 6. Electron spin echo decay of a phosphorus impurity in natural silicon (4.67% 
29 Si nuclear spins) for two different magnetic field orientations. We show experi- 
mental data from Ref. [18] together with theoretical calculations without flip-flop 
broadening (identical to Ref. [2J) and with flip-flop broadening. The latter is calcu- 
lated by directly integrating the noise spectrum shown in Fig.[5]with the appropriate 
filter function. 

are in close agreement here because for low nuclear spin density (/ = 0.0467) 
the broadenings are generally much smaller than the transition frequencies 

25-| , at least for the important pairs causing spectral diffusion. Recall that 

our theory does not account for the anisotropic hyperfine interactions. There- 
fore our theoretical results should be compared to the monotonic envelope 
enclosing the experimental data points The echo modulations due to the 
anisotropic hyperfine interaction is clearly visible at short times in the exper- 
imental data shown in Fig. [Sj These oscillations produce a loss of contrast of 
about 10% at the short time regime. Apart from this effect, the agreement 
between theory and experiment is quite good. 

Fig. [7] shows echo decay results for isotopically enriched samples (/ = 
99.23%). The experimental data is from Abe et al. [15] • Note that here the 
echo modulations are very evident, the loss of contrast reaches ~ 100%. 
The monotonic envelope on top of the experimental data is in reasonable 
agreement with the theory without broadening. However, the theory with 
broadening decays significantly faster. The difference between both theories 
increases for increasing /. This suggests that the mean-field theory proposed 
in section T4.2I overestimates the broadening. The author expects that a more 
sophisticated many-body calculation may account for this discrepancy. 



We thank Dr. A.M. Tyryshkin for pointing this out to us. 
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Fig. [8] shows the dependence of the 1/e echo decay time (T2) with the 
magnetic field angle. The shortest value of T2 is obtained when B is along 
the (111) direction (9 — 54.74°). In this case none of the nearest neighbor 
pairs have zero dipolar couplings. Only pairs i,j with parallel to the 
(100), (010), and (001) directions have their dipolar interaction turned off 
by the magic angle [fty = 54.74° implies cos 6^ = l/y/3 and fry = 0, see 
Eq. ([T0"9"J> ]. On the other hand for B || (001), T 2 is longer by a factor of 
three. This occurs because the nearest neighbor pairs, that usually give the 
strongest contribution to echo decay, are forming a magic angle with respect 
to B || (001)0 

We now discuss the time dependence of the echo envelope. The echo decay 
without broadening fits well to the expression 



<°+(2t)> = 




for a wide range of 2r centered around T2 and for all values of / (for a 
log-log plot, see Fig. 9 of Ref. [22]). Tyryshkin et al. [18] studied the time 
dependence of the natural silicon experimental data by fitting the expres- 
sion ve(2t) = e -2T / T 2e - ( 2r / T2 ) . Here T' 2 was interpreted as arising from 
a combination of spin-flip processes and the instantaneous diffusion mech- 
anism, due to the finite concentration of donors. Tyryshkin et al. reported 
= 1.1 ms and exponent n = 2.4 ± 0.1 for all sample orientation angles be- 
tween 20° < 9 < 90°. For 9 = 10° they found n = 2.6 ± 0.1, while for 9 = 0° 

11 The nearest neighbors for each site i are located at Rij = |asi(l, 1, 1), 
±asi(-l, -1, 1), ia Si (-l, 1, -1), and ±a Si (l, -1, -1). 
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Fig. 8. Orientation dependence of the 1/e echo decay time T2. 9 is the angle 
between the applied magnetic field and the crystalographic (001) direction, so that 
9 = 90° is along the (110) direction. We show experimental data for natural Si [18] 
and for isotopically enriched Si [19] . The theoretical calculations shown are without 
broadening. For natural Si, the experimental data was corrected for a 9 independent 
instantaneous diffusion decay, see Ref. [18| . Because of this the theoretical T2's 
are lower than the experimental T2's reported in Ref. [18], in contrast to what is 
observed in Fig. [6] 



n = 3.0 ± 0.2. The time dependence at angles close to the (001) direction 
is yet to be explained theoretically. At natural abundance (/ = 0.0467) the 
theory with broadening and the theory without broadening have similar time 
dependences. However, as / increases the time dependence of the broadened 
theory deviates significantly from the theory without broadening. As an ex- 
ample, for / = 0.9923 and 8 = 50° the broadened theory shows a cross-over 
from e~ T at short r < 3 [is to e~ r for r > 3 /is. This indicates that 
adding broadening to the nuclear spin excitations leads to observable effects 
in the time dependence of electron spin coherence. Unfortunately, the echo 
modulations are too strong in isotopically enriched samples (Fig. [7])- This 
makes the precise experimental determination of the time dependence of the 
echo envelope quite difficult. 

Eq. (|110[) allows us to extract scaling of the 1/e decay time T2 with the 
nuclear spin fraction /. Note that in the theory without broadening / appears 
as a pre- factor in the exponent due to p + + p_ = / 2 /4. Therefore we have 
simply 

t 2 cx r 2l2Z = r - 87 - (in) 

Abe et al. [34] measured T2 for seven isotopically engineered samples with 
/ ranging from 0.2%-100%. Their study shows that T2 must scale between 
j-0.86 anc i j-o.89 m g 00c i agreement with Eq. (|llip . 
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Fig. 9. Contribution of nuclear spins located at concentric shells around the donor 
(natural Si, 9 = 50°). Nuclear pairs closer to the center (ro < 50 A) have their 
fluctuation amplitude suppressed by the strong hyperfine field difference between 
sites i, j, forming a frozen core. The largest contribution is due to pairs located at 
50 A < r < 100 A. 



It is interesting to study the number and location of nuclear spins con- 
tributing to the noise spectrum. Fig. [9] shows the contribution due to pairs 
inside shells concentric at the donor center (for natural silicon and 6 = 50°). 
The contribution for ro < 50 A is quite small, but extends over a wide 
frequency spectrum. These nuclear spins are said to form a "frozen core" , 
because their noise amplitude is suppressed due to the strong difference in 
hyperfine fields affecting sites i, j. The frozen core of a Si:P donor has about 
3 x 10 4 nuclear spins. This frozen core effect plays an important role in other 
contexts as well such as optical spectroscopy experiments [35 . The nuclear 
spin noise theory developed here allows a quantitative description of this 
effect o From Fig. [9] it is evident that a significant fraction of the finite fre- 
quency noise power comes from the large number of nuclear spins located 
between 50 A and 100 A off the donor center (about ~ 2 x 10 5 nuclear 
spins). These pairs are satisfying a quasi-resonance condition ~ bij. 



In order to understand the frozen core effect from our analytical expression 
for the noise spectrum, assume At ^> Aj and Aij ^> bij in Eq. (|102[) . In this case 
we have cr+_ ~ Ai. From Eq. (|103[l the noise amplitude becomes ~ b\j/Ai, that is 
much smaller than bij. 
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6 Conclusions and outlook for the future 

In this chapter we described the coherent evolution of an electron spin subject 
to time dependent fluctuations along its quantization axis. We showed that 
in the Gaussian approximation the electron spin transverse magnetization 
can be expressed as a frequency integral over the magnetic noise spectrum 
multiplied by an appropriate filter function. The filter function depends on 
the particular pulse sequence used to probe spin coherence, differing substan- 
tially at low frequencies for free induction decay (FID) and Hahn echo. For a 
Gauss-Markov model (Lorentzian noise spectrum centered at zero frequency) 

_ .2 

we showed that the short time decay of the FID signal is approximately e , 
while the echo decays according to c~ r ' . 

We applied this general relationship between noise and decoherence to 
the case of a localized electron spin in isotopically engineered silicon, where 
the magnetic noise is mainly due to the dipolar fluctuation of spin- 1/2 lattice 
nuclei. The nuclear spin noise spectrum was calculated from a pair-flip-flop 
model, resulting in a linear combination of sharp transitions (delta functions). 
The echo decay due to these sharp transitions is identical to the one derived 
by the lowest order cluster expansion [21]. Next, we showed how to obtain 
a smooth noise spectrum by adding broadening to these transitions using a 
mean-field approach. The resulting noise spectrum was found to be strongly 
non-monotonic, hence qualitatively different from the usual Lorentzian spec- 
trum of a Gauss-Markov model. This structured noise spectrum is able to 

2.3 

explain the non-Markovian dynamics (~ e _r ) observed in electron spin 
echo experiments for phosphorus doped silicon. We compared the theories 
with and without broadening to two sets of experimental data, for natural 
and isotopically enriched silicon. The agreement was quite good for natural 
silicon, but not as good for 29 Si enriched samples. 

It is interesting to compare our results to the set of non-Gaussian phe- 
nomenological theories proposed a long time ago by Klauder and Anderson 
[T3]. These authors classified spectral diffusion behavior in two groups, de- 
pending on the nature of the interactions causing magnetic noise. In "Ti 
samples" the magnetic noise is caused by non-resonant spins fluctuating in- 
dividually (e.g., due to phonon emission). On the other hand the magnetic 
noise at "T 2 samples" is caused by the mutual interaction of the non-resonant 
spins (For example, a nuclear spin bath weakly coupled to the lattice is 
a "T2 sample" because the longitudinal nuclear spin relaxation time T™ is 
much longer than the transverse relaxation time T£). Klauder and Anderson 
showed that echo decay behavior in a variety of T\ samples could be described 
by a Markovian theory by making assumptions about the general shape of 
the distribution of fluctuations at any given time. While a Gauss-Markov 
model leads to echo decay of the form ~ e~ T , a Lorentz-Markov model leads 
to e~ T behavior, and intermediate non-Gaussian distributions result in e~ T 
with n between two and three. Later, Zhidomirov and Salikhov [36] showed 
that similar behavior can be obtained in T\ samples composed of a dilute dis- 
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tribution of magnetic impurities fluctuating according to a random telegraph 
noise model (Markovian with a non-Gaussian distribution). 

Nevertheless the problem of echo decay behavior in "T 2 samples" re- 
mained open. It was found empirically by many authors (see [11] and ref- 
erences therein) that echo decay behavior in "T2 samples" is usually well 

2 

fitted to the expression ~ e _r , and the Lorentz-Markov model of Klauder 
and Anderson was often invoked as a phcnomcnological explanation. Here 
we show that this behavior can be derived microscopically from a Gaussian 

model that takes into account the non-Markovian evolution of the coupled 

23 

nuclear spin bath. The resulting behavior found by us [e _r ] is not due to 
a short time approximation [the short time behavior for each pair is actu- 
ally given by e~ T , see Eq. (|96l) ]. In order to explain these experiments we 
must consider the collective finite time evolution of a large number of nuclear 
spins [note that the characteristic frequency of fluctuation for a pair flip-flop 



E-\ = 2y 6?j + A?j gets renormalized to values much larger than the dipolar 

interaction &y when the nuclear spins are subject to strong hyperfine inhomo- 
geneities = (Ai — Aj)/A}. This theoretical explanation opens the way to 
novel microscopic interpretations of a series of pulse electron spin resonance 
experiments, where the electron spin may be viewed as a spectrometer of 
low frequency magnetic noise due to a large number of nuclear spins or other 
magnetic moments. For example, a recent experiment |37j revealed that mag- 
netic noise from the surface must play a role on spin echo decay of antimony 
impurities implanted in isotopically purified silicon (with a very low density 
of 29 Si nuclear spins in the bulk). 

There are many open questions that deserve further investigation. First, 
what is the contribution of higher order nuclear spin transitions to the noise 
spectrum? This question may be answered by going beyond the simple pair 
flip-flop model assumed here, in a similar fashion as the cluster expansion 
developed in [21 , or using an alternative linked cluster expansion for the 
spin Green's function [39) . Another interesting open question is the design 
of optimal sequences for suppressing the effects of nuclear spin noise in elec- 
tron spin evolution, as was done for the random telegraph noise model in 
Ref. 40J. This is particularly important in the context of spin-based quan- 
tum computation. The efficiency of a Carr-Purcell sequence in suppressing 
the electron spin coherence decay due to a nuclear spin bath was considered 
both in the framework of a semiclassical model (see Ref. [38], where the role 
of nuclear spins greater than 1/2 was also considered) and using a cluster 
expansion approach [41 j . Recently, it was shown that the electron mediated 
inter-nuclear coupling [Eq. (|70p ] may be exploited in order to recover electron 
spin coherence lost for the nuclear spin bath [42] . We will certainly see many 
other interesting developments in the near future. 
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